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Abstract 

Digital-to-analog converters (DAC) transform signals from the abstract digital domain to 
the real analog world. In many applications, DAC's play a crucial role. 

Due to variability in the production, various errors arise that influence the performance 
of the DAC. We focus on the current errors, which describe the fluctuations in the currents 
of the various unit current elements in the DAC. A key performance measure of the DAC 
is the Integrated Non-linearity (INL), which we study in this paper. 

There are several DAC architectures. The most widely used architectures are the ther- 
mometer, the binary and the segmented architectures. We study the two extreme architec- 
tures, namely, the thermometer and the binary architectures. We assume that the current 
errors are i.i.d. normally distributed, and reformulate the INL as a functional of a Brownian 
bridge. We then proceed by investigating these functionals. For the thermometer case, the 
functional is the maximal absolute value of the Brownian bridge, which has been investi- 
gated in the literature. For the binary case, we investigate properties of the functional, 
such as its mean, variance and density. 

1 Current mismatch in digital-to-analog converters 

Digital-to- Analog converters (DAC) transform signals from the abstract digital domain to the 
real analog world. For many applications, this conversion enables the usage of the computational 
power of robust digital electronics. For example, digital audio and video, digital control, and 
telecommunications are fields that require digital-to-analog conversion. The advantageous intel- 
ligence of the applications in these fields is implemented with digital logic, e.g. microprocessors, 
and used via DA conversion in the real analog world. However, the DAC errors, at the end of 
the application chain, may decrease the performance of the whole system. Therefore, predicting 
and controlling these errors is crucial. This requirement is further emphasized in the highly 
integrated mixed-signal Systems-on-a-Chip (SoC), as we will explain now. 
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Nowadays, the SoC system solutions often integrate DAC functionality together with the 
digital logic for cost effectiveness. This requires optimal usage of the DAC resources while 
keeping the errors of the DAC within specified margins, primarily because of the two following 
reasons. Firstly, the price of the mixed-signal SoC includes the price of the DAC resources 
even when customers are not interested in the DAC functionality. Secondly, if the DAC does 
not comply with its specifications, then the entire SoC chip should be discarded. By a careful 
design, the DAC errors mainly arise from the uncertainty in the manufacturing process and, 
hence, they are random. Therefore, statistical rules are used to predict the overall performance 
for high- volume chip production. Knowledge is required that accurately links the DAC resources 
with the DAC error margins. An important example of such a relationship is the dependence 
of the Digital-to- Analog (D/A) conversion accuracy on the DAC area, i.e. the silicon area of 
the chip that is used for the DAC. Higher conversion accuracy is achieved for larger chip areas 
[21 [13]. However, too large areas introduce additional problems which may degrade the conversion 
accuracy. Thus, precise knowledge and understanding of the relationship between accuracy and 
chip area is crucial, particularly for high-volume chip systems including DAC. For a general 
introduction to DAC converters we refer to |2l[9l[18], whereas Razavi [17j also focusses technical 
aspects. 

The D/A conversion is carried out by switching certain analog quantities, such as voltages, 
currents or charges, ON or OFF. For the sake of simplicity, and without loss of generality, 
this paper will assume current as a basic analog quantity. The switching process is controlled 
by the digital input signal w G {0, 1}^, where is the length of the binary input signal. 
If all combinations of O's and I's between the digital bits produce valid input signals, i.e., 
= log2(n + 1), then the coding is called binary and is called the resolution of the DAC. A 
DAC that uses binary coding to control its current quantities is called binary DAC. The switched 
ON current quantities /„. are summed to construct the analog output signal current /out- We 
will assume that the current quantities are random and that {Iuj}]j=i are i.i.d. random variables. 
The sum of all current quantities 

n 

-^OUtmax ~ "^"J ^^■^^ 

i=i 

is associated to the maximal digital input word Wn = (1, 1, . . . , 1), and is called the full-scale 
current of the DAC. The smallest meaningful difference in the analog output, defined as the 
output least- significant bit (LSB), is defined as the full-scale current divided by the number of 
digital input codes, i.e., 

/"^OUtniax f 1 o \ 

Isb •= • U-^j 

n 

For general DAC, errors can be classified as static or dynamic. We will focus on the static errors, 
which we now introduce. 

For every digital input word Wk G {0, 1}''^, there is an analog output value /outfc- The code 
words Wk G {0, 1}^ will be assumed to be ordered. The difference between two adjacent output 
values would ideally be 7„ = ]E[JttJ, but in practice it deviates due to mismatch errors coming 
from the uncertainty of the manufacturing process, i.e. 

4ut, - 4ut,_, = 4 + (/isb - In) + DNLfc ■ Jisb, k = l,...,n. (1.3) 
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Here, Jisb — lu represents the linear error which is independent of k, and 



DNLa 



'Isb 



:i.4) 



Usb 



is the differential non-linearity in Jisb-scale. As the DAC are required to be linear devices, 
the non-hnear errors are the main concern. The integrated non-linearity (INL) measures the 
non-hnearity of the whole DAC transfer characteristic, i.e., the cumulated individual non-linear 
errors. This paper concentrates on the INL, as defined for a code Wk by 



INLfc :=^DNL, 



i=l 



-^outfc k ■ /isb 



k = 1 n, 



:i.5) 



This definition excludes the linear errors by forcing INLq = INL„ = and normalizes INL^ to 
LSB-scale. The maximal absolute INL over all codes is given by 



INL, 



max I INL A 

k=0,...,n 



:i.6) 



The statistic INL^ax is an important specification, because it indicates how linear the DAC is. 
Another important figure is YieldiNL? which indicates how many of the manufactured chips have 
INLmax under certain specified threshold. This is to say, if a DAC manufacturer guarantees 
certain linearity, then YieldiNL describes what proportion of the produced chips fall within these 
specifications. 

INL is the most popular DAC specification, see e.g. P, [18]. Its practical importance is 
very high in all DAC application fields. This general definition of INL is valid for all DAC 
architectures, and all DAC resolutions. The most commonly used architectures are the binary, 
segmented, and thermometer architectures. The segmented architecture interpolates between 
the binary and the thermometer architectures. In this paper, we consider the two extreme cases 
of thermometer and binary DAC architectures. We now describe these DAC architectures. 



For a thermometer DAC, /< 



out J. 



r(T) 

OUtfc 



is given by 



r(T) 



n. 



For a binary DAC, on the other hand, Joutt = -^out^ is given by 



N 



2™-! 



^outk ~ Bk,m luj, k — 

m=l ^=2™^! 

where the switching matrix B is an n x A^-matrix given by 

/ ■■ • \ 
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Figure 1: Thermometer DAC. The code word Wk corresponds to switching j ON ioT i < k and 
switching T^, j OFF for i > k. 
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Figure 2: Binary DAC. The code word Wk corresponds to switching those i for which Bk^. 
ON, whereas the i for which Bki = are switched OFF. The matrix B is given in (11.91) . 
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Figure 3: Segmented DAC. 
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The matrix B is constructed by writing the first n integers in (reversed order) binary coding 
into the rows of B. Note that the columns of the matrix represent the DAC bits, i.e., the 
switches of the grouped current sources. The O's represent switched OFF currents, and the I's 
represent switched ON currents. The most left column gives the switches for the least-significant 
bit current P^^l_^ = /u^, while the most right column gives the switches for the most-significant 

bit (MSB) current Iq^I^^ = Yl'j=2'^-i luy Note that 

N 2'"-l n 

Ct„ = E^«.- H 4, =5^4, =n/isb, (1.10) 

m=l j=2'"~i i=l 

SO that INL„ = 0. 

In practice, the most popular DAC architecture is the segmented one. The segmented DAC 
architectures implement one part of the input digital bits in a binary way and the other part 
in a thermometer way. That is how the advantages of both the binary and the thermometer 
architectures can be combined. 

Figure [3] shows a segmented DAC architecture. The LSB part is implemented in a binary way 
and the MSB part in a thermometer way. The output current for a given level of segmentation 
S is expressed by 

JV-S 2™-l [A:/2^-'5j (m+l)2^-'5-l 

m=l j=2"'^i m=l j=m2^~S 

The parameter S* G {0, . . . , A^} determines the interpolation between the binary and the ther- 
mometer parts of the DAC architecture. Indeed, for = 0, the segmented DAC architecture 
is transformed to a fully binary architecture. On the other hand, for S = N, the segmented 
DAC architecture is transformed to a fully thermometer architecture. Note that in both extreme 
cases the same unit current sources are used. Thus, the performance difference is only due 
to the way the J„. are combined to construct the output current /outfe- In the binary DAC, 
the unit currents are first grouped and then switched ON or OFF, while in the thermometer 
DAC, the unit currents are individually switched ON or OFF. More detailed discussion on DAC 
architectures can be found in the hterature [9], HH 119] . 

As discussed, due to manufacturing related mismatch, the currents lu^ always deviate from 
their designed values, i.e., /„. = + ej. The mean value is chosen by the DAC designer, 
whereas Ej is the random error due to the manufacturing process which we model as an i.i.d. 
sequence of normal random variable with zero mean and variance cr^. Nevertheless, our results 
remain valid when the {ej} are i.i.d. distributed with sufficiently many moments. The ratio 
(Ju/Iu is known as the relative current matching, i.e., the unit current matching. The relative 
matching determines the required transistor area for the particular manufacturing process. The 
smaller Oujlui the more accurate l^p but the larger the required area. For more details on the 
dependence between relative matching and transistor area we refer to [H]. 

Once (Jujlu is specified, the required transistor area can be calculated. However, the rela- 
tionship between Oujlu and DAC non-linearity, which is crucial to determine the proportion of 
chips complying to specifications, has never been determined analytically. So far, DAC engineers 
have used either Monte Carlo simulations or approximations. 



5 



The Monte Carlo simulations of a DAC model produce empirical results to suggest some 
design specifications Oujlu- Although not accurate, this approach is very practical, see [5]. A 
problem arises for the design of high-resolution DAC, for which N is large. The complexity of 
the DAC model increases by a factor of two for every additional bit in resolution. For example, 
for = 14, ra = 16383 unit elements have to be simulated. Therefore, Monte Carlo simulations 
are not practical for higher resolutions, because they become complex and slow. 

On the other hand, a number of analytical approximations can be found in the literature. The 
analytical attempts to describe INL, and in particular INLmax? started in 1986. The approach 
in [11] disregards the correlation between the DAC outputs for different input codes. Bastos [2] 
proposes a much simpler formula, which considers only the deviation of the transfer characteristic 
at the mid-scale DAC output, which can be a rough, though simple, estimation of the INLmax- 
Another approximation was given by van den Bosch et al. [1] by assuming that, if the DAC 
static transfer characteristic at any code INL has error equal to the target value, e.g., INL^ = 
O.SJisb, then there should be a 50% chance that ultimately INLmax is smaller than the target 
value, i.e., INL^, < 0.5 /isb for all k. The major approximation inaccuracy is in the probability 
that both the positive and negative INL limits are reached for the same DAC sample, e.g., 
INLfc < — 0.5/isb is disregarded. Though this approach derives a convenient normal distribution 
for INLjnax; it is inaccurate for higher resolutions, as we show in more detail in this paper. In 
general, approximations lead to transistor over-design, i.e., a transistor area that is too large. 
In this paper, on the other hand, we shall present an exact analytical formula, for which no 
approximation is necessary. 

Due to the lack of exact analytical formulation of INL and the high complexity of DAC 
model simulations, the statistics used in industry for a high volume chip production are hard to 
predict. Here we think of the statistics YieldiNL, the INLmax distribution, INL^ax deviation and 
mean. Furthermore, the advantages of some redundancy based approaches relying on the DAC 
statistical INL properties cannot be theoretically estimated, see e.g. [TB]. Finally, up to now 
the main DAC architectures, i.e. binary, thermometer, and segmented, cannot be distinguished 
with respect to their static linearity properties, so they are wrongly considered identical [2[ H] . 
One conclusion from our results is that the INL for binary and thermometer architectures are 
different 

Implications of the results derived in this paper for the field of DAC's can be found in a 
second paper A comparison with the results in [U [5l [11] is summarized in [151 Table 1]. 

2 Thermometer coding: Maximum of a Brownian bridge 

For the thermometer coding, we can describe the INL as functional of a Brownian bridge as 
follows. 

Theorem 2.1 (INLmax for the thermometer coding). As n —>■ oo, 



and L"^, where the limit X is characterized by 




(2.1) 



in distribution, in 



X = max \Bt 

te[o,i] 



(2.2) 
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for a Brownian bridge {Bs}s&[o,i\, o-nd 



E[X] = In 2 0.86875, Var(X) = — - ^(In 2)^ 0.0677732, (2.3) 

and 

oo 

P(X <x) = l + 2^(-l)^e-2'='"', x>0. (2.4) 

Note that in (12. ip . we multiply by /isb rather than by /„. For the convergence in distribution, 
this makes no difference what so ever, since /isb converges to a.s. by the strong law of large 
numbers. However, the convergence in and fails if we multiply by /„, since in this case, the 
expected value of INL^ = (/out^ — k ■ Iish)/I\sh is not defined, whereas |INLfc| has infinite mean. 

We recall that a Brownian bridge {Bs}s£[o,i] is a Markov process on [0, 1] that is obtained 
from a Wiener process or Brownian motion in either of the following two ways: 

(Bl) 

B, = Ws-sW,, 
where {H/s}sg[o,i] is a Wiener process; 

(B2) 

Bs = Wi, 

where {VF^'jsgp,!] is a Wiener process conditioned on Wl = 0. 

For an introduction to Wiener processes and Brownian bridges we refer to [7]. For the equivalence 
of (Bl) and (B2) cf. P pp. 83-85]. 

Proposition 2.2. For the thermometer coding, 

max \f3 -kIisb\=(Tuy/n max \Bk/n\ , (2.5) 

k=0,...,n ' *" ' k=l,...,n ' ' 

'D 

where {Bt}telo,i] is a Brownian Bridge process and = denotes equality in distribution. 

Proof. Since {/« }j=i,...,n is a family of i.i.d. normally distributed random variables with mean 
/„ and variance cx^, we have that 

k 

/isb = (4,-/^ 4) (2.6) 

is normally distributed with mean and variance ka"^. Hence for a Brownian motion {VFt}t>o, 

J<S, - fc/isb = (T^Wk = a^^/^Wk/n, (2.7) 

where we used Brownian scaling in the last distributional equality. By substituting k = n and 

— _____ 
using (11.21) we obtain n(/isb — nlu) = Ou^nW\. Combined with (12. 7p . this yields 

/outfc - k /isb = cTuVn I Wk/n - -Wl j = auVnBk/n (2.8) 



for a Brownian bridge process {Bs}s&[o,i], where we used (Bl) above. After taking absolute 
values on both sides of (12.81) and the maximum over k = 1, . . . ,n, we obtain (12.51) . □ 

The maximum over the discrete time points k/n, k = 1, ... ,n, in (12.51) can be replaced by 
the supremum over the whole interval [0, 1] by using the following lemma: 

Lemma 2.3. For C > A, 



P 



max -Bfc/n — max IBA 
k=l,...,n ' ' te[o,i] 



> C./^ 1 < 4ni-^'/« + 2n-^^"/«. (2.9) 



n 



In particular, maxfc=i^... „ |-Bfc/n| converges to maxjgp^i] \Bt\ in distribution as n ^ oo. Moreover, 
L}- and L"^ -convergence holds. 

The proof of Lemma 12.31 is deferred to the appendix. We now use Lemma 12.31 to complete 
the proof of Theorem 12.11 

Proof of Theorem \2.1[ The convergence in (12. ip follows from Lemmas 12.21 and 12. 3[ For the 
probability of the upper tail of X (12.41) we refer to [3, (11.39)], the computation of mean and 
variance of X is straightforward integration. For example, we have that 



E[X] 



/ F{X>x)dx = 2 V(-l)'~'e-2'='^'dx = 2V(-l)^-i^ = ^ln2. 



(2.10) 

The interchange of summation and integral can be justified by looking at = X V e. □ 



3 Binary coding: Block increments of a Brownian bridge 
3.1 Results and overview proof 

In this section we prove the following theorem characterizing the binary coding statistic. 
Theorem 3.1 (INLmax for the binary coding). As n ^ oo, 

-^INL^a. — . M, (3.1) 
in distribution, in and L^, where the limit M is characterized by 



M = - 5^ 2-('+i)/2 



2 

1=1 



(3.2) 



for a family {Zi}'^^ of i.i.d. standard normal random variables. The expectation of M is given 
by 

oo 

E[M] = {27i)-^/^ {2-^ - 2-^^) ^''^ ^ 0.839792, (3.3) 
1=1 

and the variance Var(M) is computed explicitly in (I3.23P below and can be approximated as 

Var(M) ^ 0.08007. (3.4) 
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Note that in (13. ip . we multiply by Jisb rather than by /„. As explained for Theorem 12 .11 this 
makes no difference for the convergence in distribution, though the convergence in and 
fails. 

The proof of Theorem 13.11 is organized as follows. In Lemma 13.31 we prove the convergence 
in (13. ip . where the limit is characterized in terms of increments of Brownian bridges. After this, 
Proposition 13.41 shows that the weak limit M can be expressed as the weighted sum of standard 
normal variables, which proves (13. 2p . Lemma [3.71 states the expression for mean and variance of 
M. Finally, we give an approximation to the density of M in Section 13.51 



3.2 A Brownian bridge representation of the binary INL^ax 

Our aim is to derive an expression for INLmax for the binary coding. First we express the 
non-linearity of the current steering DAC in terms of a functional of a Brownian bridge. 



Lemma 3.2. 



N 



k=l,...,n 

where {Bs}1^q is a Brownian bridge. 

Proof. Let {Wt}tG[o,oo) be a Wiener process, then 

V 



(3.5) 



m=l 



- Iu = OTuVn {Wj/n - W(j-i)/n) 



j = l,...,n. 



(3.6) 



V 



V 



where = represents equality in distribution. We will further write = rather than =, because we 
are interested in the distribution only. Furthermore, 



1 " - 1 1 

hsh - ^« = - (^«, - In) = - (^vWu = (^uWi 
n ^-^ ^ ' n \ n 



(3.7) 



where we recall that the are i.i.d. Miju-, o'„)-distributed. 

We want to calculate IomI^, ~ ^h&h for k being a power of 2 first. The advantage is that, for 
k = 2*^-^ only the mth block {J„,|j = 2^"-^ . . . , 2™ - 1} contributes: 



t(B) nm-l T 

-'out„^_i ^ -'Isb 



2™-l 



/ / 2"" — 1 

= auVn ( W v^-i Wi 

\\ " n 

= au\/n ( i? 2'"-i — i? 2"i-i-i 

V " n 

where, in the last line, we used the representation (Bl) of Brownian bridges. 



2m-i _ I 

2m-i_i Wl 

n 
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For calculating maxfe=i,...,„ |/outi ~ ^-^isb|, we need to take the maximum over every configu- 



ration of contributing blocks, hence 



max 

k=l,...,n 



-'outfc fi'-'lsb 



Ou^Jn max 
ac{i,...,Af} 



mG3 



(3. 



:=A/jv 



We denote by 3i* the subset of {1, . . . , A^} for which the maximum in fl3.8p is achieved, and use 
the abbreviation 



m = l,...,N, 



(3.9) 



for the increment of the Brownian bridge on the interval [ 



2m-i_i 2™-!" 
n ' 



'-^] ■ Without loss of gener- 
ality, we may assume that X^mea* Bn,m is positive, otherwise the same argument for —B holds. 
Clearly, 

me 3* ^ Bnm>0. 



Furthermore, 



N 



N 



N 



- i?i - i?o - ^ Bn,m - ^ Bn,m '^Bn,m>0 + X] ^K^^ 



,<o- 



m=l 



m=l 



m=l 



Using fl3.10p in the first line and (13.111) in the second, we obtain 



N 



M, 



N 



m=l 

I ^ . 1 - . 

m=l 

9 E P"'^ 



iV 



m=l 



!<0 



m=l 



which is fl3.5D. 



We write 



as in (13. Sp . and define 



M 



N 



max 



auy/n k=l,...,n 



Isb 



^ oo 

M:= 152-0-1) 

where 5 denotes a Brownian bridge. 

Lemma 3.3. There exists a constant C > 0, such that 



P 



Mn - M 



>e]< 



2-7V/2 



(3.10) 



(3.11) 



□ 



(3.12) 



(3.13) 



for every e > 0. In particular, M^r converges to M in distribution as N 
converges in and in . 



(3.14) 

oo. Moreover, it 
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We show in Proposition 13.41 below that M = M. The proof of Lemma 13.31 is deferred to the 
appendix. The combination of Lemmas 13.21 and 13.31 yields the convergence 

oo 

-5^152-. -52-,+i)| (3.15) 
1=0 

in distribution, in and L"^, as N —>■ oo (and thus also n = 2^ — 1 ^ oo). This proves fl3.ll) . 



max 

n k=l n 



-kh 



Isb 



3.3 A representation of M in terms of i.i.d. standard normals 

In this section, we will prove the following representation formula, which expresses M in terms 
of independent standard normal random variables. 

Proposition 3.4 (Rewrite of M in terms of standard normals). Let Zi, Z2, ... he a sequence of 
i.i.d. standard normal distributed random variables. Then, M can be expressed as 



^ 00 



2 

1=1 



(3.16) 



In other words, M in (13.131) has the same distribution as M in (13.21) . 

In order to obtain the limit law of M, we have made use of the representation (Bl) of 
the Brownian bridge. Now we will primarily use (B2). We will essentially use the following 
well-known property of Brownian motion: 

Lemma 3.5 (The conditional law of the middle point). Let {Ws}%,q be a standard Brownian 
motion. Then, the distribution of conditional on Wt = z is a normal distribution with mean 
I and variance |. 

Lemma 3.6 (The distribution of {-B2-'}zSi)- The distribution of {B2-i}'^i is equal to 

I 

^2-. = 5^2-(^'+i)/2-('-^')Z„ (3.17) 
i=i 

where {Zj}'^^^ are independent and identically distributed standard normal random variables. 

Proof. The proof is by induction in /. For Z = 1, we use Lemma [3.51 together with the fact that 
{BsW^Q is a Brownian motion conditioned on Bi = 0. This implies that the distribution of Bi 

is equal to a normal random variable with mean and variance |. Therefore, 

= \z,, (3.18) 

where Zi is a standard normal distribution. This initializes the induction hypothesis. To advance 
it, assume that 

i?2-a-i) = Yl 2-(^'+i)/2-(/-i-J-)^^.. (3.19) 
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Then, again using Lemma I3.5[ the distribution of conditionally on i?2-(i-i) is a nor- 

mal distribution with mean l/2 52-('-i) ^ind variance 2^^'+^^ Therefore, if we denote Zi = 
2('+i)/2(^2j2-! — |i?2-(!-i)), then Zi is a standard normal random variable independent of i?2-{i-i)- 
As a consequence, we have that 

1 ' 

B,-^ = 2-(^+^)/2^; + -52-a-i) = J2 2-(^'+i)/2-a-,)^^.^ (3^20) 

i=i 

where in the last step, we have used the induction hypothesis. □ 



Proof of Proposition \3.4\ As a consequence of Lemma 13.61 we obtain the identity 

i-i 

52-a-i) - 52-. = J2 2-(^^+i)/2-a-,)^^. _ 2-('+i)/2z, (3.21) 

Thus, by f l3:T3D . we obtain §J6^. □ 
3.4 The moments of M 

In this section, we identify the first two moments of M: 

Lemma 3.7 (Moments of M). The expectation of M is given by (13.31) . that is 



^ CO 

[M] = —= {2-^ - 2^2') ^ 0.839792, (3.22) 



E 



1=1 

and the variance by 



Var(M) = IJ2 (2^^ " 2"^^/^ (S"^' - 2"^^)^/^ ( ^/w^ - I + Pik arcsin(p.,)) + ^1 " ^) 

(3.23) 



i<«<fc 

with 



_2-{l+k) 

plk = , (3.24) 
The variance of M can be approximated by Var(M) ^ 0.080066. 

In the proof of Lemma [37fl we will make use of the following property of the bivariate normal 
distribution: 

Lemma 3.8 (The expected product of absolute values of normals). Let Yi, Y2 be two standard 
normal random variables having a bivariate normal joint distribution with correlation coefficient 
p. Then 

E[|ri||r2|] = -v^r^P^+ — arctan ( / \ = - f ^/l - + p arcsin(p)') . (3.25) 

TT 71 ViyX_p2/ 7]- V / 
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For a proof of Lemma [3.81 see e.g. flUi Excercise 15.6]. 

Proof of Lemma \3. 7\ The expression for the mean is easily derived. We note that 

i-i 

Ni := B^-o-i) - = 2-(^'+i)/2-(/-i)^^. _ 2-a+i)/2^, Z = 1, 2, . . . (3.26) 

i=i 

is normally distributed with mean and variance 

I 

VI = 2-(^'+i)-2('-^') = 2-' - 2-2', (3.27) 
and rewrite M in ([32D as M = i YT=i \^iV 

Since, for a normal random variable Z with mean 

and variance cx^, we have that 

m\] = \J^. (3.28) 
the representation formula (13.161) allows us to identify the mean of the random variable M as 

E[M] = - ^ E [liVd] = ^ E = ^ E (2"' - 2-^0 ^ 0-839792. (3.29) 
For the variance of M we expand 

Var(M) = ^(2 ^ Cov(|iV,|, |iVfc|) + ^ Var (|iV,|) j . (3.30) 

\ l<l<k 1=1 / 

The variance term is not too hard, as 

Var (|iV,|) = (E[iVf] - nm?) = " ^) ^^.31) 

and therefore, 

f; Var (m) = (1 - 1) f;(2- - 2-'") = 5 (1 - f ) . (3.32) 

1=1 ^ ^ 1=1 ^ ^ 

Fix now 1 < Z < A;. Then 

Cov(iV,, Nk) = Gov ( E 2-(^'+^)/2-(z-i)^^. _ 2-(^+^)/^Zi, ^ 2-(^'+^)/2-(z-i)^. 
l-l 

= E 2-(^'+^)-('--'')-(''--'') - 2-('+i)-('=-') = _2-('+'=). (3.33) 
i=i 

Therefore, {Ni,Nk) is a bivariate normal distribution with mean (0,0), variances (vi,Vk), and 
correlation coefficient 

Cov(iV,,iV,) -2-('+^) 
v/Var(Arz)v/Var(Arfc) ^2^' - 2-2' V2-*=-2-2'=' ^ ' ^ 
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Denoting (1^, Yk) a bivariate normal distribution with means 0, variances 1 and correlation 
coefficient pik, we have 

Cov(|iV,|, \Nk\) = V^Cov(|l^|, mi) (3.35) 
= y^{E[\Yi\\Yk\]-E[\Yi\]E[\Yk\]) 

= ^^(E[\Yi\\Yk\]-^y 
By Lemma [3.8[ as well as fl3.30p - fl3.35p . we obtain 

Var(M) = IY1 [\l^-pl -^+Pik arcsin(pifc)) + 1 M - ^ j . (3.36) 

l<l<k ^ ^ 

Using fl3.27l) and fl3.34l) we can approximate Var(M) numerically, which yields fl3.4p . □ 
Having completed the proofs of fl3.ip - fl3.4l) . the proof of Theorem 12.11 is complete. 



3.5 Approximating the density of M 

We now derive a formula for the density of M. Without loss of generality, we may assume that 
= and (Ttt = 1, i.e., the are standard normal distributed. 

By denoting N = (iVi, N2, . . . ) and Z = (Zi, Z2, . . . ), we rewrite (13.261) with the help of the 
(infinite) matrix L as N = —L ■ Z, where 

.2-a+i)/2-(/-i) if J < i- 

L,, = <j 2-('+i)/2 iij = l- (3.37) 

if j>/. 

Note that L is a lower triangular matrix. 

We will approximate the density of the infinite sum M = 1/2 Yll^i ^"7 the finite sum 



M^-") = -^|iV;|, mGM. (3.38) 



2 

1=1 

Writing N*"' = (A^i, . . . , Nm), Z'™' = (Zi, . . . , Z^), and L'™' for the upper left m x m corner 
of the infinite matrix L, we have that N'™' = — L*™' ■ Z^""\ In particular, N'™' is normally 
distributed with mean (0, . . . , 0) and covariance matrix S'™' = L'™' ■ (L'™')^, where we write 
^^{mj-jT fQj, ^YiQ transpose of the matrix L'""'. Note that 

(S""> ., = < , „ 3.39 

Given the mean and covariance matrix of a multivariate normal distribution, its density is 
known to be 



14 



where n = (ni, . . . , Um) G M™. We write IN'^'I for the pointwise absolute value {\Ni\, . . . , \Nm\) 
of the m-dimensional vector N'""'. Its density is given by 

Wl(n)= E 7^4^ .det E(-))V2 {-^(^n)(S''"))-^ {anf\ , n G [0, 00)-, 

(3.41) 

where (cm) = {(Titii, . . . ,0-^12^)- See e.g., [U (6.3.20)]. 

The determinant det S'"' and the inverse (S*"')"^ of the covariance matrix are easy to com- 
pute since L*'"' is a triangular matrix. The results are stated in the following two lemmas. 

Lemma 3.9. For all m G N, the determinant ofJ]^""^ is 

detS'™' = 2-™(""+3)/2. (3.42) 

Proof. First we note that 

detS*™' = (detL""')l (3.43) 

Since L^™^ is a triangular matrix, its determinant is obtained by multiplying the entries on the 
diagonal, and hence 

det S'™' = ( J] -L<7^) = ( J] 2-('+i)/2) = 2-™("^+3)/2. (3.44) 
i=i j=i 

□ 



Lemma 3.10. For all m G N, i/ie inverse ofE^""^ is given by 

, , ,,„i f(j A/)2(-'+')/2-i for J I; 

hi ={r^Lj-i i (3-45) 

Proof. Since L'™' is a triangular matrix, it is easy to see that the inverse (L''"')~^ is given by 

r 2(^-i)/2 if J > I- 

(L""')-! = ! 2(^+i)/2 if J = /; (3.46) 

[0 ifj</. 

For j < I we have that 

iA« j 
(Se-));/ = ^^((Le"')-!),, ((L<-')-i),, = J2 2(^--i)/22('-i)/2 = J 2(^-+')/2-\ (3.47) 
fc=i fc=i 



whereas for j = I, 



(S(™')7/ = 5^((L<'"')"')?fc = 5^ 2(^-i)/2 + 2^j+m = + 1) 2^-1. (3.48) 

fc=l A:=l 

□ 
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In order to calculate the density of M'™' = |(|A''i| + ■ ■ ■ + |A^m|) at ?/ > 0, we have to integrate 
(13.401) over the (m — l)-dimensional surface {(ni, . . . , rim) I 2?/ = |^i| + ■ — h |^m|}- This leads to 
the following formula: 



2y 2y-ni 2y-rn nm-2 

fMim)iy)=2j j ■■■ j /|N(m)|(ni,n2, . . . (2y - ni n„_i)) 



drim-i ■ ■ ■ dn2 drii, y G [0, oo). 
Note that there are m — 1 integrals. Thus for m = 1 there are no integrals and 



(3.49) 



fMwiv) = 2/|N(i)|(2y) =4^/-exp{-8y"}, y e [0,oo). 



(3.50) 



Finally, (I3.49P gives us a formula for the density of M*'"^ The only (numerical) problem are 
the integrals. For larger values of m, an intelligent way of numerical integration seems necessary. 



2.5 



1.5 



0.5 




0.25 0.5 0.75 1 1.25 1.5 1.75 

Figure 4: The density of M'""' for different values of m, in contrast to the density of X. 

However, for sufficiently small m, a mathematical standard package, such as Mathematica, gives 
good approximations (see Figure HI). 



3.6 A disintegration approach to the density of M 

In this subsection we present a different approach to the density of M. We define the quantity 



^ oo 



Wo 



(3.51) 



1=1 
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for a Wiener process {W,},e[o,i]. It is obvious that M = M\^^^^, cf. flCTD and (B2) on page [71 
Let / be the Fourier transform of the joint distribution of M and Wi, i.e. 

f{ki,h) =E[exp{iikiM + k2Wi)}]. (3.52) 

Using the independence and stationarity of the increments of the Wiener process in fl3.5ip we 
obtain 

oo 

Hh, k^) = TTE[exp{2(A;i |W^2-a-i) - W^2-'l + k2 (H^2-o-i) - W^-i))}], (3.53) 

T=i' ' 

=hi(ki,k2) 

where 

hi{ki,k2) =E[eicY>{i{ki\Z\ + k2Z)]\, (3.54) 

and Z is a A/'(0, 2~')-distributed normal random variable. Once we have computed hi, we obtain 
the joint density of (M, Wi) via inverse Fourier transformation as 

fM,Wji^^y) = J j expj-ifcix - ik2y} /(fci, ^2) '^^^^^l ■ (3-55) 

The density of M equals the density of M conditioned on Wi = 0, whence 

fiviix) = ^ = V27[ [ [ exp{-ikix} f{ki, ^2) ^f^^ ^ . (3.56) 



fwM J J ^ ^ .-^^j,^., ^2vr2)2 

It remains to calculate hi{ki,k2). If we let 

hi{k)=E[e''^l{z>o}]^ (3-57) 

then hi{ki, ^2) = hi{ki + ^2) + hi{ki — k2). The dependence of hi on / can easily be eliminated 
by scaling: 

hi{k) = h{2-^/^k), (3.58) 
where h is as in fl3.57p with Z being A/'(0, l)-distributed. It can be shown that 

h{k) = -e-''"'^+i-^ e^'/2^x, (3.59) 



2 2v^ 
but the inverse Fourier transform in (I3.56p seems intractable. 



4 Conclusions 

We have derived the distribution of the INLmax in terms of Brownian bridges. This distributional 
identity holds for all architectures, in particular also for the segmented one. 

In the thermometer case we have identified the limiting distribution of INLmax as the absolute 
maximum of Brownian bridges, which is well-known in the literature. 

For the binary case, we have identified the limiting distribution of INLmax in terms of a 
Brownian bridge. We have further provided a representation in terms of independent standard 
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normal variables and have computed the mean and variance of INLmax- Finally, we have given 
a procedure that approximates the density. 

We want to emphasize that the INLmax in the thermometer case and in the binary case behave 
differently. Although the densities look alike, e.g. the upper tails are quite close to each other, 
there are significant changes in the lower tail. The thermometer case has, compared to the binary 
case, a slightly larger mean, but a slightly smaller variance. Even though the distributions in 
the two cases are close, they are not the same. 

We still miss the distribution function for the binary case. Random sums of the type in fl3.13p 
have appeared in the literature. In particular the quantity 

oo 

S = Y,^-'Vi (4.1) 

1=1 

where {V;}^]^ are i.i.d. exponential random variables arise in a variety of applications, see e.g. 
Ott, Kemperman and Mathis [13], Guillemin, Robert and Zwart [8J and Litvak and van Zwet 
|12j . The density of S can be expressed in terms of an infinite sum, cf. ^L3\ Section 5]. When 
the summands are independent uniform random variables, i.e., {V;};^^ are i.i.d. uniform random 
variables on (0, 1), the density of S can be computed explicitly [B]. 

Furthermore, it would be interesting to extend the results to the segmented case. In partic- 
ular, it would be of interest to investigate which limiting INL^ax distribution has the smallest 
mean. This should correspond to the optimal DAC architecture. Practical implications of our 
results can be found in a companion paper [15]. 



A Proof of Lemmas 12.31 and 13.3 

Proof of Lemma \2. 3[ We first prove (12. 9p . Let {Bs}seYd,i] be a Brownian bridge. Then, 



max -Bfc/n — max \BA < max \Bk/n — BA 

k=l,...,n ' ' ie[0,l] k=l,...,n ' ' 



(A.i: 



Using representation (Bl) above, we can further bound ( lA.ll) from above by 

max \Wk/n-Wt\ + -\Wi\ 

k=l,...,n ' ' n 



(A.2) 



for a Wiener process {VTsjsgp,!]. Using the Markov property and Brownian scaling, we obtain 

(A.3) 



that for A; = 1, . . . , n. 



max \Wk/n - Wt\ = —= max \Wt\, 
\hizl ii ' ' ^Jn te[o,i] 
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where = stands for equality in distribution. Hence, for C > 0, 



max -Bfc/n — max IBA 
k=l,...,n ' ' te[o,i] 



>^ /logn 



n 



<PI ma..; ma. |».,,„ - »'.| > ^;./!2Sii I + P Vll > 



k=i,...,n I tgfiiiii ' 2 \ n \ \ n 2 \ n 

< nP^^max iVTtl > ^^/lognj + f(^\Wi\ > ^a/ti log nj . (A.4) 



(A.5) 



For the first term, we bound for every 6 > 0, 

P( max \Wt\ >b) < 2P( max Wt > b] = A¥(Wi > b) < 4e"^'/^ 
Vte[o,i]' J V*e[o,i] J 

where we use the reflection principle [TJ Theorem (6), p. 526] in the second and a standard bound 
on the tail of standard normals in the third step. Substituting b = C/2 y/logn we obtain 

nP^^max \ Wt\ > ^v^logn^ < n4n-(^/2)V2 ^ 4^i-c^V8_ (^.6) 

For the second term in (lA.4p . we obtain analogously 

^(^\Wi\ > ^v^n logn^ < 2f(wi > ^^nlogn^ < 2^"^'"/^ (A.7) 

The bound f[01) . together with (jAj])-(jAjj) proves ([23]). 
For the convergence in L^, we use that, with 

Xn = max l-Bfc/nl , X = max \BA , (A. 8) 
k=l,...,n ' ' ' te[o,i] 

we have that 

E\Xr,-X\= F{\Xn- X\>t)dt. (A.9) 



We split between t < 4./^ and t > 4./^. For t < 4. r^, we bound P(|X„ - X| > t) < 1 



while for t > 4^^, we use (lAlej) and (IATTI) to bound 

P(|X„ - X| > t) < 6ne-"*'/l (A.IO) 
Substitution of these bounds yields 

E|X„-X|= rF(\Xn-X\>t)dt<Axf^+6n H e'^'" dt = O { xl^^) . (A.ll) 

V n 

The convergence in follows similarly, now using 

E[{Xn-Xy]= 2tF{\Xn- X\>t)dt. (A.12) 
Jo 

We leave the details to the reader. □ 
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Proof of Lemma \3. 3[ We observe that, by (13.51) . fl3.12p and flS.lSp . and replacing m hj N — m, 



/n|+ E 



2-m 



N-l 

m=l in=N 

Therefore, for an arbitrary chosen constant e > 0, we have that 



(A.13) 



P 



Mm - M 



\m=l 



Bo 



\m=N 



^ e 



{II) 



(A. 14) 

Using representation (Bl), we see that for a Brownian bridge {i?s}sg[o,i], and any < s < t < 1, 
the following holds for every constant C > 0: 



¥{\Bt-B,\>C)<¥{\Wt-Ws-{t-s)W,\>C) 

<^{\Wt-Ws\ > C/2) + F{\{t-s)Wi\ > C/2). 



Using the Markov inequality and Gaussian scaling we obtain 



4 2 



P (15* - 5,1 > C) < ^Vt^s E[\W,\] + -it - s) E[\W,\] <-^=-{t~s) 



(A.15) 



(A. 16) 



We use ( ]A.16I) to bound (/) in (1A.14I) from above by 



m=l 



B, 



(2N-m -l)/n 



> 



2(iV- 1) 



m=l 



\{N-l) 



< 



£2^/2' 



e V 2™ 2^-1 

(A.17) 

which converges to as ^ oo. The second term (//) in (IA.14P can be bounded using the 
Markov inequality by 



(//) < - E 



^ ] |i?2-™ — B2 



m=N 



(A.18) 



Using (I3.29p . this expectation can be computed as 



E 



^00 " ^ 00 

— / LBo-m — Bn-m-A = — -= y (2 
2 ^ ' ^ 



m=N 



(A.19) 



m=N 



and 



(//) <'-^±(2- 



r < 



m=N 



e^{V2-r 



2-7V/2 



(A.20) 



which converges to as A^ ^ cxo. 

The combination of ([04]) . (IATTTI) and ([OOl) shows that, for C = 8 + 4(^27? - 



P 



Mn - M 



> e] < 2-^/2 



(A.21) 



for every e > 0, i.e., M^r converges to M in distribution as A^ — > 00. 

The convergence in and in follows as in the proof of Lemma 12. 3[ 



□ 
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